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Abstract 

t*^ ' The chiral Abelian Higgs model contains an interesting class of solitons found by Rubakov and 

o ■ 

CN ' Tavkhelidze. These objects carry non-zero fermion number Np (or Chern-Simons number Ncs^ 

^ , what is the same because of the chiral anomaly) and are stable for sufficiently large Np- In this 



paper we study the properties of these anomalous solitons. We find that their energy- versus- 



' fermion-number ratio is given by ~ ^cs o^" ^ ~ ^cs depending on the structure of the scalar 

(N 

^ ■ potential. For the former case we demonstrate that there is a lower bound on the soliton energy, 

which reads E > c n"^^ , where c is some parameter expressed through the masses and coupling 
\ constants of the theory. We construct the anomalous solitons numerically accounting both for 



Higgs and gauge dynamics and show that they are not spherically symmetric. The thin wall 
approximation valid for macroscopic solutions with Nqs ^ 1 is discussed as well. 
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I. INTRODUCTION 



Solitons - stable localized solutions to the classical equations of motion of non-linear field 
theory represent an interesting class of particle-like states in quantum field theory. Well 



HQ, 



known examples include topological solitons such as ;he iink in 1 + 1 dimensions 
the vortex in 2 + 1 dimensions 0, 4|, the monopole s], [g] and the skyrmion [7] in 3 + 1 
dimensions. The stability of these solutions is ensured by topological reasons. Another class 
of solutions - non-topological solitons or Q-ba lls j^, 9], are stable because of conservation 
of some global Abelian [9i] or non-Abelian 10| charge. The above mentioned objects exist 
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in pure bosonic theories. Yet another type of solitons, "electroweak bags'' 
fact that fermions can be trapped inside a spherical cavity created by a (otherwise unstable) 
space-dependent scalar field. 

In 13|, referred to as RT in the following, a new class of solitons was found, for which 



the presence of chiral fermions and Abelian gauge symmetry play the key role. In very 
general terms, the construction is based on the following observation. Consider an Abelian 
Higgs model with chiral fermions and arrange the Yukawa couplings in such a way that all 
fermions get a mass rriF from the Higgs condensate. Then, the energy of Np well separated 
fermions is Ep ~ mpNp- Due to the chiral anomaly these fermions can be converted to a 
gauge configuration carrying the Chern-Simons number Nqs = Np, and the energy of the 
gauge-Higgs system E with non-zero Nqs may appear to be smaller than Ep. If true, then 
one gets a stable soliton characterized by Ncs- Indeed it was demonstrated in RT using a 
variational principle that the energy E of the system with Ncs 7^ is bounded from above by 
EqN^^q, where Eq is some constant. Thus the bosonic configuration is always energetically 
more favorable than the collection of Np fermions at sufficiently large Ncs- The stable 
bosonic configuration will be referred to as "anomalous Abelian soliton" in the following. In 
'act, anomalous Abelian solitons have some common features with Hopf solitons, studied in 

The Abelian character of the gauge group is essential for absolute stability of anomalous 
solitons jisl . 1^. Indeed, if the gauge group is non-Abelian, fermions can always be converted 



into a gauge vacuum configuration 



17, hi 



[19I . which may carry arbitrary integer Chern- 
Simons number. In other words, non-Abelian anomalous solitons, if they exist, can only be 
metastable. 
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Interestingly, anomalous solitons may potentially exist 15|, ll6| in the Standard Model, 
since it contains an Abelian U(l) gauge group and chiral fermions. This issue, however, has 
not been clarified yet. 

This paper is devoted to the detailed study of anomalous Abelian solitons. We demon- 
strate, with the use of different inequalities from functional analysis, that for a certain class 
of scalar potentials also a lower hound on the soliton mass, E > c N^^g, can be established, 
where c is some parameter related to the masses and coupling constants of the underly- 
ing theory. We construct the anomalous solitons numerically, taking into account both the 
dynamics of the gauge and the Higgs field, and find that the solitons are not spherically sym- 
metric. We also discuss the thin wall approximation valid in the limit of large Chern-Simons 
number, which will allow us to remove the Higgs dynamics from consideration. 

The organisation of the article is as follows. In the next section, we review the Tavkhelidze 
a.d Rubakov co.st..uct,o,> Q Q of anomalous .oUto... o.de. to „,ake the paper self- 
contained. In Sec. IIIII we derive a lower bound on the energy of the solitons and discuss 
their general properties. In Sec. [IV] we present the numerical solutions of the field equations 
as a function of Chern-Simons number and coupling constants of the theory. Section |V] is 
devoted to the solution in the thin wall approximation. Finally, we summarize our main 
results in Sec. |VT1 



II. THE RUBAKOV- TAVKHELIDZE SOLITON 



A. The model 



Consider an Abelian Higgs model with the complex scalar field $, a pair of left-handed 
fermions with opposite charges ip^ and a pair of neutral right-handed fermions ip}^ and 
V^^, and chiral interaction between the fermions and the U{1) gauge field. The Lagrangian 
is 

-Ai [i^ii^],^ + h.c.) - A2 {^li^l^* + h.c.) (1) 
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with the potential 

V{\^\) = X{\^\'' -v^ (2) 

and Dfj^ = 9^ — igA^. The vacuum expectation value (VEV) of the scalar field is equal 
to V and the gauge and Higgs bosons obtain the masses my = \plgv and mu = 
respectively. For simplicity we choose the two Yukawa couplings to be equal to each other, 
Xp = Xi = X2, such that all fermions have equal masses mp = Xpv. Note that we are forced 
to introduce at least two fermions in order to make the theory free from gauge anomalies. 
The total fermionic current 

2 

i=l 

is anomalous in this model, 

where / = 2 is the number of left-handed fermions. 



B. Instability of fermionic matter 

In this subsection we follow fisl. [isl in order to explain why the fermionic matter becomes 
unstable at sufficiently high fermionic density. 

At zero fermionic density, the ground state of the boson fields is A = and I'I'I = v, which 
is called the normal state, following RT. Consider now a homogeneous system of infinite 
volume containing fermions with number density np. The system is supposed to be neutral 
with respect to the gauge charge. We only consider the weak-coupling limit (A ^ 1 and 
g"^ <^ 1), where we can treat the fields and $ as classical condensates. The neutrality of 
the system implies Aq = and that the fields A and $ are time- independent. The fermions 
can be characterized by the chemical potential fip, which is, at zero temperature, the energy 
up to which the Fermi levels are filled. The fermion number density np is related to the 
Fermi energy fip hj 

nF = f 



2 



377 

We will use the unitary gauge Im($) = and denote = Re($) in the following. The 
fermions can be integrated out 20|] and, up to corrections of order jjp^, the static energy 
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functional of the bosonic fields is given by 



£x , (4) 



Ncs[A] = J-^ / e'^'F,,Ak d'x = I A ■ B d'x (5) 



where B = 'V x A is the magnetic field. The first four terms in (j4]) are the classical energy 
density of bosons, while the last (Chern-Simons) term is due to the interaction with fermions. 
The physical reason for the appearance of the Chern-Simons term is as follows. As the gauge 
field A increases, some fermionic energy levels cross the zero-energy line and the number of 
real fermions decreases by fNcs, where 

-,2 /" 2 

327r2 J 167r2 

is the Chern-Simons number of the gauge field. 

The quadratic part of the static energy functional (jlj) has a negative mode for 

fflp > /icrit , 

where 

/icrit = :r^V ■ (6) 

Thus the normal ground state of fermionic matter is absolutely unstable at 

/^crit 



n, > = ^ . (7) 



The negative mode is given by 

e(a:;) = (ei cos(fc • a;) + 62 sin(A; ■ a;)) , (8) 

where k = my and ei.2 are real polarization vectors orthogonal to each other and to k. Note 
that a perturbation A{x) = ae{x), where a is a small amplitude has the Chern-Simons 
density 
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Ties = nA ■ B = -ka^ . 



As the amplitude of the unstable mode grows, the term in the energy acts as a 

positive mass term for the Higgs field, which leads to the disappearance of the Higgs field 
condensate. The system undergoes a transition to a state with = containing a negligible 
number of real fermions and a gauge field condensate A ^ with non-zero Chern-Simons 
number. We will call this state the abnormal state, as in RT. 



C. Domain of abnormal matter 



Still following 13|, [iSj, we show in this subsection, that a finite domain of the abnormal 
state with sufficiently large Chern-Simons number has a lower energy than a system contain- 
ing Np = fNcs real fermions. The stable configuration of gauge and scalar field condensate 
is the minimum of the static energy functional for the bosonic fields 



E\A, 



2\2 



(Tx 



under the constraint of a constant Chern-Simons number ([5]). Varying the functional 

2 

E[A,ct>]-^^-^Ncs[A] 



(9) 



(10) 



with respect to A and where /i is a Lagrangian multiplier, we get the following field 
equations 

V^^ - /A2(/) - 2A(02 - t;2)0 = , (11a) 
V X (V X A) -2/i(V X A) + 2/0^^ = . (lib) 

Solutions of Eqs. (ITT]) with finite energy have to satisfy the boundary conditions 



(j) V and A ^ for ^ 

For a static solution {Aci{x),(f)ci{x)} of Eqs. (ITTl) we obtain 

167r2 



oo 



E[A^ 



ch Vd\ 



/i- 



-N, 



9' 



cs 



X 



(12) 



where the first term in Eq. iQ has been integrated by parts and flllbp has been used. 

Let us consider a compact domain of finite volume V in the abnormal state, embedded 
into the normal vacuum. We suppose that the size R ~ V^^^ of the domain is much larger 
than the two characteristic length scales in the model, the thickness of the domain wall 
~ l/mn and the penetration length ~ 1/my of the magnetic field. Then, the contribution 
of the scalar field to the energy (the second term of Eq. f|T2l) ) can be written 



(13) 



where ly = V/R^ is the shape factor for the domain V. The surface energy, which is of 
order R'^v^, can be neglected. Since we want to minimize the total energy. 



1 2 

E = fx^Ncs + Xv^R' Iv , 



(14) 



with respect to R, we have to find how the Lagrangian multipher /z scales with R. Since 
= inside the domain of the abnormal state, it follows from Eq. flllbl) that fi ~ 1/R. 
Thus we can write /i = v/{2R), where u is a. numerical coefficient. The total energy as a 
function of R becomes 

E{R) = Svr^^^ + Xv'R^Iy (15) 

g R 

and minimizing E[R) with respect to R leads to 

\ 3Iv } mv\ l3 } ' ^ ' 

E^V-2[^y'(IvPrjNi',\ (17) 

where we expressed R and E in terms of the vector boson mass my and the parameter (3 
defined by 

,3 = 4 = ^. (18) 
Both the size of the domain R and the total energy E grow slowly with increasing Ncs- If 
Ncs > NcTit, the energy (fTTIl is smaller than the energy of the normal state with Np = fNcs 
real fermions for which ~ nipNp. Parametrically, A^^cnt ~ (p^)^- Therefore, the domain 
of abnormal matter is stable. Since R ~ v^^N^g, we get R ^ v^^ for large Ncs, and the 
surface energy R'^v^ ~ "^^cs becomes indeed negligible at large Ncs-, as we assumed above. 



In [l5| the domain was considered to be a sphere of radius R and Eq. (lllbp was solved 
analytically (remember that = inside the domain). The solution can be expressed in 
terms of Bessel functions and one gets 

where is the first node of the Bessel function J3/2. Inserting z/ = ^0 and Jy = 47r/3 into 



Eqs. ([TSD and ([HD we find exactly Eqs. (6.7) and (6.8) of [15|. 

This completes our review of the main results in RT, who proved the existence of anoma- 
lous Abelian solitons. The remainder of the present paper is devoted to the detailed study of 
the structure of these solitons. In particular, we would like to answer the following questions: 

• What is the structure of the gauge and Higgs field condensates of anomalous solitons? 

• What is their exact shape? Are they spherically symmetric? 
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Is the exponent 3/4 in the power-law dependence of the energy on Ncs universal for 
anomalous solitons? 



• Does the structure of the solitons depend on the choice of the scalar potential? 

• What happens at small Nqs, when the dynamics of the scalar field has to be taken 
into account? 

III. PROPERTIES OF THE SOLITON 
A. Lower bound on the energy 

In the following we derive a lower bound on the energy of anomalous solitons, which reads 

E > cN'J,' , (19) 

where c is a constant depending on the VEV of the Higgs field and the coupling constants. 
Such a bound on the energy is rather unusual because of the fractional power of the topo- 



logical charge. A bound of this type was first obtained by Vakulenko and Kapitanksy 21 1 
for solitons of the non-linear a-model, for which the topological charge is the Hopf number. 
The Vakulenko bound was later improved (i.e. increasing the constant c) by Kundu and 



Rybakov |22| and Ward [23 1. 



The derivation of ( fT9l) is based on the use of some inequalities of functional analysis, 
which we review in appendix El Let us start with an estimate of the Chern-Simons number 

(ISD, 



Ncs = g'ci j A-Bd'x 

< g^c, (^j d^x^ (^j |A|6 ci^^y^' , (20) 

where Ci = l/(167r^) and we have used the Holder inequality (]A2|) for p = 6/5 and g = 6. 
The first factor on the r.h.s. of fl20|) is estimated using flA3p . 



5/6 / ^ \^/^ / /• \^/^ 

\Bf^d?x] <[ \B\^d?x] [ \B\(fx] , (21) 



and the second factor with the use of the Gagliardo-Nirenberg-Sobolev inequality (lASp 

1/6 / /• \ 1/2 



(Tx) < C2 / V|A| (Px . (22) 



Rosen 2J] found the smallest possible value for the constant C2 in the previous inequality 
to be 

The integrand of (1221) is further bounded from above, 

3 3 

\A~ 



2/3 

C2 = -^(-] 



i=l *>i=i 

3 

= |V X + ^ diAjdjAi . (23) 

Imposing the Coulomb gauge V ■ A = 0, the last term in Eq. fl23p vanishes through integra- 
tion by parts and we obtain 

j \V\A\\^ d^x < j \B\^ d^x . (24) 
Inserting ^ and ^ into ^ yields 

Ncs < 9^ciC2 / d^'x] / |B| d^'x , (25) 



or equivalently, 

N^s < 1^1' d^x^ (^j \B\ c^^x^ , (26) 

where 

C = {ciC2f'^ = ^ . (27) 

At this point, it is worth mentioning [23], where it was claimed that the constant C can 
even be reduced to C = 1/ (2567r^) . If we define an average magnetic field by 

we can write down an inequality for the magnetic energy 

12 1 / D\ 1/2 / N 3/4 



(29) 



In order to obtain an inequality for the total energy rather than for -Emagn, we consider 
the scale transformation 

0A(a;) = 0c/ (Aa;) , 
Aa{x) = AA,i{Ax) , 
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where {(j)ci{x), Aci{x)} is a static solution of the field equations with finite energy. This 
transformation leaves the Chern-Simons number unchanged. The static energy 

E{A) = E[Aa{x),Mx)] 
must have a minimum at A = 1. We get for the energy 



where y = Ax. Requiring that 

dE 
dA 

leads to the relation 



= 

A=l 



j^Sy = j [(V0,,)' + g'Alct^l + 3^(0,0] d'y , (30) 

which is exact and has to be satisfied by any static solution of the field equations. From 
( |30|) we deduce 

E = ^i5;^agn + IJ [(V0eO' + /^^^ d'x (31) 

and using (l29l) we get 



which is a lower bound for the energy of anomalous solitons, provided the average magnetic 
field B is bounded from below. We will derive such a bound in the next subsection using 
known results from Ginzburg-Landau theory of superconductors. 



B. Anomalous solitons and superconductors 



The static energy (Q is just the relativistic version of the Ginzburg-Landau free energy 

= V corresponds to the supercon- 



for a superconductor (see e.g. [2g]). The state with 
ducting state, while = coincides with the normal-conducting state. It is well known 
from superconductivity theory that the superconducting state is destroyed by large external 
magnetic fields and that the way in which this superconductivity breaking occurs, depends 
on the ratio /3 of the Higgs mass to the gauge boson mass squared (cf. Eq. (|T8l) ). 

In superconductors of the first kind ((3 < 1) the superconducting state can persist only 
up to a critical magnetic field (see e.g. 



26[), which is given by = {2Xy/^v^. At B > B^ an 
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intermediate state is formed, in which normal-conducting domains (with = and B > Be) 
emerge within the superconducting state {(p = v, B = 0). Thus for type I superconductors 
the lower bound on the average magnetic field is provided by the critical magnetic field, i.e. 



B>B,= ^ml 



and we finally get the lower bound 



E>^C-^l^i5^l^'^NTs for (5<1. (33) 

On the other hand, for superconductors of the second kind (/? > 1) the magnetic field pen- 
etrates the superconductor along vortex filaments with quantized magnetic flux (Abrikosov 



vortices 



3)- 



For B(,\ < B < Bc2 the superconducting sample is in a mixed state forming 
a lattice of Abrikosov vortices. The lower critical field Bd marks the point, when the first 
vortex filaments appears and is given by 

$0 

where e is the energy per unit length of one vortex filament and $o is the unit flux quantum. 
For (3^1 the vortex energy can be calculated analytically and the lower critical field 



becomes 



25| 



^ci = Alnff). (34) 



2v^ V2, 

If we suppose, that for /? ^ 1 the state with the lowest energy is a lattice of vortices of total 
length L, we can write the average magnetic field as 

- f\B\^d^x 2Le ^„ 
j \B\ d?x L$o 



and we obtain for the lower bound 



3 



C. Shape of the solitons 



-/2 r //QMi/2 
E > 1±C-^I^ 



ln(^ 



2 



'^nT, for (36) 



From (macroscopic) superconductivity theory we can also gain an insight on the shape 
of anomalous solitons at large A^cs- Consider a surface separating a normal- conducting 
and a superconducting domain within a type I superconductor in the intermediate state. 
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Because in the superconducting domain B = 0, the continuity of the orthogonal component 
B± across the surface tells us, that the magnetic field in the normal-conducting domain 
has to be tangential on this surface. Furthermore the magnetic field has to be equal to the 
critical field B^. Therefore we expect, at least for /3 < 1, the magnetic field to be tangential 
and equal to on the boundary of the solitons. However, it is not possible to contruct 
a constant tangent vector field on a sphere. More generally, on a closed surface of genus 
there is no continuous tangent vector field, which is non-zero everywhere on that surface 
(Poincare-Hopf theorem). The simplest manifold on which such a vector field can exist, is 
a torus (the surface has to be of genus 1). Consequently, the highest possible symmetry of 
anomalous solitons is axial symmetry. 



D. Dependence on V{(l)) 

Up to this point we have discussed anomalous solitons for the potential (I2j). In this 
subsection we consider solitons for different scalar potentials. We shall see, that the 3/4 
power-law of the soliton energy, E ~ A^,^^ , is not universal, but depends on the scalar 
potential. 

Let us consider three different classes of potentials shown in Fig. [TJ For potentials with 




FIG. 1: Three choices for the scalar potential V{(\)): (a) Potential with V{i^ = 0) = Vq > (non- 
zero critical magnetic field), (b) Potential with V{j^ = 0) = (zero critical magnetic field), (c) 
Potential with V{(^ = 0) < (the ground state = u is metastable). 



1/(0 = 0) = Vo > 0, such as the potential ([2]) (see Fig. 1(a)), there is a non-zero critical 
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magnetic field, which is given by 



2 



l^n >0 . 



(37) 



The critical magnetic field is the energy density of the scalar field in the abnormal state 
= 0. On the boundary of the two phases, (j) = v and = 0, there is a pressure balance 
between the pressure of the magnetic field and the vacuum energy density. The surface 
contribution to the energy of the scalar field is negligible. 



On the other hand for a potential with Vq = (Fig. 1(b)), the critical magnetic field 
vanishes. In this case the energy of the Higgs field is dominated by a surface term and is 
given by 



where I a is a numerical coefficient. Instead of Eq. ( |T5|) we obtain 

E{R) 



57r 



9' R 

Minimization with respect to R shows that 

/ \ 1/3 
1/6 9i^\ 



(38) 



(39) 



R 

E 



3(327r8)V6 



/J 

A 



-N, 



1/3 



my 



CS ' 



T ,.2\ 1/3 
1 



my 2/3 

2 ^^CS 



(40) 
(41) 



9 J 9 

In the next section, we will numerically verify this result for a specific potential of the type 



of Fig. 1(b) 



Finally, for a potential with Vq < (see Fig. 1(c) ), there is no critical magnetic field. The 
normal ground state = w is metastable and anomalous solitons do not exist in this case. 



IV. STRUCTURE OF ANOMALOUS SOLITONS 



In this section we construct anomalous solitons numerically. We present numerical solu- 
tions to the field equations as a function of the Chern-Simons number Ncs and the parameter 
(3. The structure of the Higgs and gauge field condensates is analyzed in detail. Finally we 
discuss the dependence of the solutions on the scalar potential. 
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A. Numerical Construction 



In numerical calculations it is convenient to work in dimensionless coordinates x = myx, 
and with dimensionless fields (p = (p/v and A = A/v. In terms of these quantities the field 
equations (ITTi) are 



Vx (VxA)-2/2VxA + ^A 








(42a) 
(42b) 



where /x is the (dimensionless) Lagrangian multiplier. We want to solve Eqs. f l42l) subject to 
the constraint of constant Chern-Simons number, 

1 



Ncs\A] 



327r2 



A-B d^x 



(43) 



which determines the Lagrangian multiplier /i. 

It followed from the discussion in section IIII CI that the highest possible symmetry for the 
solitons is axial symmetry. Therefore we restrict ourselves to axially symmetric solutions 
and use cylindrical coordinates (r, z), for which Eqs. fH2|) are reduced to a system of 
non-linear PDEs in the rz-plane: 



- dr (r dr<p) + dl<P - 

r 


'\{Al + Al + A-:)<P-^-{^-- 


1)0 


= 0, 


(44a) 




- dr {rA^) 
r 


d^ {d^Ar - drA^) - 2/i d^A^ - 


fAr 


= 0, 


(44b) 




+ dlA^ + 2/i {d,Ar - drA,) - ^ 


h'A^ 


= 0, 


(44c) 




^ dr [r {drA, - d-Ar)\ + ^ {dr {t A^)\ - 


fA^ 


= . 


(44d) 



For axial symmetry we have 



^cs = ^ y '^^v i^^Ar - drA^) drdz 



(45) 



after integration by parts and over the angle if. Note that we are omitting hats from now 
on. 

Solutions of Eqs. (jHj) are calculated by discretizing the fields 0(r, z) and A(r, z) on a 
rectangular box [0, L] x [—Z^ Z] in the r2;-plane and then solving the corresponding system 
of non-linear equations using Newton's algorithm. The boundary conditions on the edges of 
the box are 



1 



and 



1 



(46a) 
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for the Higgs field and 



and 



r=L 



z=-Z 



z=Z 



(46b) 



for the gauge field. On the axis of symmetry we imposed the regularity conditions 



dr 



drA, 



r=0 







and Ar 



r=0 



A,. 



r=0 



. 



r=0 



(46c) 



For a detailed description of the numerical procedure we used to solve the PDEs, we refer 
the reader to appendix [Bl 

B. Dependence of the soliton energy on Ncs and (3 

We calculated solutions for 1 < Ncs < 1-2 ■ 10^ and 0.4 < /5 < 16. In Fig. [2] we show the 




50 100 



FIG. 2: Logarithmic plot of the soliton energy as a function of Ncs foi' f3 = 1- The dots show the 
energies calculated from the numerical solutions for Ncs = 30 to Ncs = 12'000. The solid line is 
a fit of the form E = clNq^ with a ~ 118.826 my / . The dashed line corresponds to the lower 
bound (1331). 



energies obtained from the numerical solutions as a function of Ncs for = 1- The energy 
satisfies E = aN^g , with a ~ 118.83 my j ^ which confirms the 3/4 power-law dependence 
on Ncs- The numerical solution gives an energy about three times larger than the lower 
bound (I33D, which is ^ 39.7 my j g'^ N^J^ fo^ = 1- 

Fig. [3] shows the soliton energy as function of (3 at Ncs = 1'600. We find E 3.01 • 
10"^ my j with p 0.23. The (5^1^ dependence is expected from Eq. ( IT7|) (or from (l33l) ) 
for /5 < 1. 
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FIG. 3: Logarithmic plot of the energy as a function of (3. The dots represent the energies of the 
numerical solutions for /? between 0.4 and 16 (for Nqs = 1'600). The solid line is the fit -E = bpP 
giving p ~ 0.23 and b ~ 3.007 • 10^ my /g^ . The dashed line corresponds to the lower bound (j33p . 

All the solutions we found consist of a single domain with <C w and \B\ > sur- 
rounded by normal vacuum. For small Chern-Simons number Ncs < 20 the soliton has a 
size of the order of the thickness of the domain wall (~ l/mn), but becomes larger when 
Ncs is increased. 

This is the behavior we expect for (3 < 1 (cf. intermediate state of type I supercon- 
ductors). For (3 > 1 one could imagine solutions containing a network of closed vortices 
similar to Abrikosov vortices in type II superconductors, which could have a lower energy 
than our solutions for /3 > 1. Therefore we point out, that our solutions for (3 > 1 might 
be metastable with respect to the decay into a vortex-type solution of equal Chern-Simons 
number. However, we did not find any solutions of this type. 

C. Structure of the Higgs and the gauge fields 

Let us now look at the scalar and gauge field configurations of the solitons. We will 
concentrate on the specific example with (3=1 and Nqs = 12'000. 

Figs, m and [5] show a 3D-plot of the scalar field in the r2;-plane and a 2D-plot of 
along the axis of symmetry r = 0, respectively. The thickness of the domain wall (~ 1/m//), 
which separates the two phases -C f and = f , is small compared to the size of the soliton 
(thin wall approximation). 

Moreover Fig.lHshows that the soliton is not spherically symmetric. We can determine the 
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FIG. 4: The scalar field in the rz-plane for Ncs = 12'000 and /? = 1. Note that (j) is shown in 
units of V and the coordinates r and z are in units of l/my- The rectangular box in the rz-plane 
used to obtain the solution in this case was [0, 52] x [—40,40]. 

1 [ . 



oL — ...■■■...■■■I-...--..-..-.; 

-40 -20 20 40 

z 

FIG. 5: The scalar field 4> along the axis of symmetry r = (again for Nqs = 12'000 and /? = 1). 
The dots show the values of (p on the grid points of the numerical solution. The resolution is 
0.8/mv. 

shape of the soliton numerically. For this end, we define the contour of constant scalar field 
= 0.5 f as the boundary of the soliton and fit a 2-dimensional surface to that boundary. 
This is done in Fig. El in which we show the contour = 0.5 v in the rz-plane. The contour 
follows very precisely a circle of radius ~ 27.132/m^/ centered at (IS.lGS/my, 0). The 
contour deviates from a perfect circle just in a narrow region around the axis of symmetry, 
whose size is of the order of a few l/rriv- Only in this small region, the fields behave 
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FIG. 6: The dotted line shows the contour of constant scalar field (p = 0.5 in the rz-plane. The 
solid line is the best fitting circle to the dotted line, which is the circle around (18.163,0) with 
radius 27.132 (in units of l/my). 

microscopically. Motivated by the amazing numerical result of Fig. El we conjecture that 
asymptotically for large Ncsi the soliton has the shape of a so-called spindle torus. A spindle 
torus is the surface of revolution obtained by rotating a circular arc around an axis, where 
the arc is intersecting that axis. However, we were not able to prove this result analytically. 
We now turn to the gauge field. As depicted in Fig. [TJ the magnetic field is confined to 



|B| 




FIG. 7: 3D-plot of the magnetic field in the rz-plane in units of the critical magnetic field Be- 
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the region where the scalar field is vanishing or small. Inside this domain the amplitude of 
the magnetic field is larger than the critical magnetic field Be- 

On the boundary \B\ = and B is tangential to the domain boundary (see Fig. [H]). 
Remember that this behavior is expected from (type I) superconductors in the intermediate 

1 

0.8 ; ' 

B| I . 6 - ; 
0.4 - • 

0.2 ; ; 

10 20 30 40 

r 

FIG. 8: The dots show the tangential component of the magnetic field along the boundary 
of the domain as a function of the distance r from the axis of symmetry (units for Sy: Be = 1). 
Sufficiently far away from the axis of symmetry the magnetic field is constant and equal to the 
critical field Be, which is represented by the solid line. 

state: on the boundaries of normal-conducting regions the magnetic field is tangential and 
equal to B^. Outside this boundary, where the scalar field starts to deviate from zero, the 
magnetic field decays exponentially, because in the broken phase the gauge boson becomes 
massive. 

To get an idea of how the magnetic field lines behave, we show in Fig. [9] a plot of the 
"toroidal" component B^ and in Fig. [10] the "poloidal" component Bp = {Br, B^) of the 
magnetic field. The toroidal component of the magnetic field is maximal at (20.8/my,0) 
and vanishes on the boundary of the spindle torus, where the field is completely poloidal. 
The poloidal field is maximal at the origin and zero at (27.4/my,0). As it was already 
shown in Fig. [8|, \Bp\ = B^ along the boundary of the spindle torus. In short terms, the 
poloidal component of the magnetic field wraps around the toroidal component, giving rise 
to the non-zero helicity of the magnetic field (which is the same as Chern-Simons number 
for Abelian fields). 
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FIG. 9: Contour plot of the toroidal component B^p of the magnetic field (in units of the critical 
magnetic field Be). The shown contours are from to 2.4 with a step size of 0.4 between the 
contours. The maximal value of i?^ is located at (r, z) ~ (20.8, 0) and is equal to 2.4 Be- 

Finally we show in Fig. [TT]tlie energy density of the soliton [Ncs = 12'000 and (3 = 1). 
The different parts of the energy are 

^magn = J ^ d'x ^ 1.00158 ■ 10^ my / 9^ , (47) 
Ed = J [(V0)2 + ^^AV^] d^x ^ 3.985 ■ 10^ my/g'^ , (48) 
Ev = j A(0^ - v'^f d^x ^ 3.2123 • 10^ my / g"^ ■ (49) 

The virial theorem fl30l) is very well satisfied (-Emagn = 3i?v + to an accuracy of less 
than 1%). The part E-q of the energy is a surface term and is negligible for sufficiently large 
Chern-Simons number. 



D. Solitons for different scalar potentials 



In section UlI Dl we have seen that for a scalar potential like the one shown in Fig. 1(b) 
the energy of anomalous solitons is ~ A''^'^. 
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FIG. 10: The poloidal component Bp of the magnetic field. The sohd hnes are projections of the 
field lines onto the rz-plane. The colors indicate the modulus of Bp (in units of Be). The poloidal 
field is maximal at the origin and vanishes at (27.4,0). On the boundary of the spindle torus 
\Bp\ = Be- 
in order to verify this conjecture, we also solved the field equations for the potential 

F(0)=7(0-^^)V. (50) 

If we again use dimensionless fields and coordinates as before, the field equations fj42|) for 
the potential (jSOll read 

- 1 - f [0 (0 - 1)' + 0' (0 - 1)] = , (51a) 

V X (V X A) - 2/iV X A + (p^A = , (51b) 

with 13 = 2'~f/g'^. The constraint equation (H3|) remains unchanged. We solved Eqs. (15T!) for 
10 < Ncs < 700 and (3 = 1. In Fig. [12] we show the energy of the obtained solutions as a 
function of Ncs- We find 

E ^ 94.866 mv/g^ N^g , (52) 

which confirms the claim of section IIII D[ Fig. [13] shows the energy density the soliton with 
Ncs = 500 and (3 = 1- The main contribution still comes from the magnetic field, but one 
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r 40 

-40 

FIG. 11: 3D plot of the energy density for the soliton with Ncs = 12'000 and /3 = 1 (in units of 
the vacuum energy density Xv^). The maximum of the energy density is at the origin and is equal 
to £max ~ 12.56 Xv^ = 6.28 
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FIG. 12: Logarithmic plot of the energy as a function of Nqs for anomalous solitons obtained with 
the potential (jSOp . The dots show the values of our numerical solutions and the solid line is the 
best fit of the form E = aN^^^, which gives a ~ 94.866 my/g'^. 

can clearly see that there is a contribution located on the surface, which is the energy of the 
scalar field. As well in this case, the solitons are not spherical. 

V. THIN WALL APPROXIMATION 

In this section we construct the anomalous solitons in the thin wall approximation valid 
for large Ncs- In this limit the solitons are described by a compact domain V, inside which 
= 0. Outside the domain we have (j) = v and A = 0. In order to determine A inside V we 
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FIG. 13: The energy density of the sohton with Ncs = 500 and /3 = 1 in units of 7^^/16 (which 
is the potential barrier between the two minima of the potential). 



have to solve Eq. (Illbp for = 



V X (V X A) = 2/i (V X A) 



(53) 



where A = on the boundary dV. In principle, the surface dV has also to be determined 
by minimization. But because we have found numerically in the previous section, that the 
solitons have the shape of a spindle torus for large Ncs, we take dV to be a spindle torus 
of fixed size described by two parameters Rg and > Rs (see Fig. [H]). We solve Eq. (|53|) 




FIG. 14: Cross-section of a spindle torus in the rz-plane with radii Rg and r^. The spindle torus 
is obtained by rotating the circle centered at {Rs,0) with radius > Kg around the z axis. 



inside the spindle torus as a function of Rg and r^. Then we minimize the soliton energy 
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with respect to these two parameters. We use dimensionless cyhndrical coordinates {g, ip, () 
defined by 



X = RsQCOSif , 

y = RsQsmip , 
z = RsC. 

In these coordinates (!53l) simphfies to a system of hnear partial difi^erential equations for the 
potentials Ag, A^ and in the gC^-^Aane 



{d^A, - d,A 



-dlA^ 



-de {qA^) 



-de[g{dcA,~dM] 
Q 



V {d^A, - d^A^) , 



V -dg {qA^) 
Q 



with V = 2fiRs. 

Integration of (I54al) and (I54cp gives 



(54a) 
(54b) 

(54c) 



d^A, - OA. 



RgB^p 



and using (155!) in (I54bp leads to 



do 



Q 



. 



The Chern-Simons number becomes 

Ncs 



9' 



A-B (fx 



167r2 

2 

2 / AipBip d X 



9 



Al (fx 



(55) 



(56) 



(57) 



M'^Rs 

where the second line is found by partial integration and Eq. fl55|) has been used on the 
third line. We solved Eq. 0561) numerically with = on the boundary of the spindle 
torus and under the constraint flFTl) . Because Eq. fISBl) is a linear eigenvalue problem, the 
solution which minimizes the energy is the eigenfunction to the lowest positive eigenvalue 
i/^, normalized with Eq. (!57|) . It turns out, that the eigenvalue i/ is a function of the ratio 
Qs = fs/Rs, which is shown in Fig. [151 The total energy of the solution is given by Eq. (fT5|) 
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FIG. 15: The eigenvalue as a function of the parameter Qg. 

for R = Rg. After minimization with respect to Rg, the size of the spindle torus Rg and the 
energy are given by (cf. Eqs. (fT6l) and ( fTTl) ) 

1/4 



-Rule's 









SIviQg) 


my 





IQTT^ulgg 



3/4 



2 ^'CS ' 



where the shape factor for the spindle torus is 



Air 



Ivies) = tt'qI + 2ny^^^ + y (^'^ - 1) + 2^^^' 



arcsm | — 



(58) 
(59) 



(60) 



In Fig. [TB]we show the function E{gg), which can be minimized numerically with respect to 
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FIG. 16: The total energy as a function of the parameter Qg. The minimal energy is at Qg ~ 1.495 

^3/4 



and is ^min ~ 119.06 P'^/^{mv / g^)Nl^^ 
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Qs- There is a minimum at Qs ~ 1.495, for which we get 



^ 1.742 



1 fN, 



1/4 



1/4 "^V ^.3/4 

9^ 



E ^ 119.065 (3 —^^cs 



(61) 
(62) 



These results are in excellent agreement with the results of the previous section, where we 
obtained 



R, ^ 1.726 



1 fN, 



niv V P 



ws 



1/4 



E ^ 118.826 /3V4!p:^rV4 _ 



(63) 
(64) 



The macroscopic solutions discussed in the present section coincide with the solutions of the 
previous section for large Ncs- To illustrate this, we compare in Fig. [17] the magnetic fields 



|B| 




FIG. 17: The amplitude of the magnetic field \B\ in the plane z = of the soliton with Nqs = 
12'000 and [5 = 1. The dotted line shows the solution of the previous section, the solid line is the 
solution in the thin wall (macroscopic) approximation discussed in the present section. The units 
are Be = 1. 

of the macroscopic solution and the solution of the previous section. 



VI. CONCLUSIONS 



In this article we studied anomalous Abelian solitons, a class of solitons of the chiral 
Abelian Higgs model found by Rubakov and Tavkhelidze. These solitons carry a non-zero 
Chern-Simons number Ncs and they are stable for large Ncs- We showed that their energy- 
versus-fermion-number ratio is given by ~ ^c^s or ~ -^c's depending on the structure 
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of the scalar potential. For the former case we derived a lower bound for the soliton energy, 
E > cNJg , where c is a constant depending on the coupling constants and the VEV of the 
Higgs field. 

We discussed the structure of anomalous solitons by solving numerically the field equa- 
tions for both the gauge and the Higgs field. It turned out that the solitons are not spherically 
symmetric. For large Chern-Simons number they have the shape of a spindle torus. In this 
limit we also constructed the solutions in the thin wall (macroscopic) approximation. 

At the moment it is not clear, if (metastable) anomalous solitons exist in the standard 
model. It will be interesting to extend the present study to the full electroweak theory with 
SU{2) X U{1) gauge group. 
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APPENDIX A: INEQUALITIES 

In this appendix we review the inequalities from functional analysis, which have been 
used in Sec. IIII Al to derive the lower bound on the energy of the anomalous solitons. 

1. Holder inequality 

Let be an open subset of Mf" and p, g > 1 such that 



Then, for all functions / G U{VL) and g G L'^iP), the product fg G V-{VL) satisfies the 
Holder inequality. 



1 



1 



+ 



1 . 



P 



\\fg\W < VIlAMl. , 



(Al) 



where the norm on U'{VL) is defined by 
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Note that for p = q = 2 the Holder inequahty reduces to the Cauchy-Schwarz inequahty. 
For vector fields u{x) : Q ^ M.'^ the L^-norm is defined by 

i/p 

I'^W LP = \ I (u, uY^"^ dx 



where (■, ■) denotes the scalar product of vectors in M°'. Thus for vector fields the Holder 
inequahty reads 



< WuWt.pWv 



LP W^^WLi 



(A2) 



2. "Modified" Holder inequality 



As a consequence of flAip we have for s > 1 that 



i= < ll/llir|l/llir^ ^feu{n). 



The proof of flASP is very simple: defining the functions g and h by 

h=\fr-\ 

we get 



(A3) 



11/11^ = \\9h\l 



2~s 



\9 



1/(2-^) dx 



s-1 



dx 



2-s 



2a-2 
2 



I/I dx 



\f?dx 



lA\f\\%-' , 



where ( lAlD has been used on the second line. 
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3. Gagliardo-Nirenberg-Sobolev inequality 



The Sobolev spaces W'^''^{Q), where p > 1 and n G N, are the normed spaces of functions 
defined by 

W'^ifl) = {/ e LP{Q)\ V« G |a| < n : 9^/ G ^^(fi)} . 
For 1 < p < d we define the Sobolev conjugate of p by 

dp 

P = 1 • 

a — p 

Then, there is a constant c{d,p), such that 



<c{d,p)\\Df\\L. , yfew''''{n), (A4) 

where D is the weak derivative. (1A4I) is known as Gaghardo-Nirenberg-Sobolev (GNS) in- 
equahty. For f2 = R'' the best possible constant in (lA4p is 

i/d 



c{d,p) 



p> 1 . 



-1 


d—p 


1/p 


r(d+i) 


l/cZ 


-p 


d{p-l) 




r(d/p)r(<i+i-d/p)|S'*-i 





Here, |5"^ ^| denotes the surface area of the {d — l)-dimensional unit sphere. 



r{d/2) 



T{d) • 

In Sec. nil Al we have used (lA4p for continuosly differentiable functions with compact support 
and {d,p) = (3,2). In this case the weak derivative is just the ordinary gradient and the 
GNS inequality reads 

1 / 2 \ 

i6<c(3,2) ||V/|U2 = — M IIV/IU^, WfeCliR'). (A5) 



APPENDIX B: NUMERICAL PROCEDURE 

In this appendix we explain the numerical procedure which was used to obtain numerical 
solutions of the PDEs (l44l) subject to the constraint (H5l) and boundary conditions (H6l) . 
Solutions are calculated on the rectangular box [0, L] x [—Z, Z] in the rz-plane. The dimen- 
sions of the box L and Z have to be at least a few times N^^^, since the expected size of the 
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soliton is i? ~ ^cs- box is divided into a grid of A^^, x A^^ points (rj, Zi) given by 



Ti = hr{i — 1) , for i = 1, . . . , N, 

h 
~2 



z^ = -l{2j-l-N,) , for j = l,...,N,, 



where hr = L/{Nr — 1) and hz = 2Z/{Nz — 1). Because the coordinates {r,z) are in units 
of I/my and we expect the fields to vary on scales of order I/uih and and 
have to be smaller than (3^/"^ to ensure a sufficient resolution of the solutions. Therefore, the 
number of grid points has to be increased going either to larger Ncs oi to lower (3. As an 
example for Ncs = 12'000 and /5 = 1, we were using the rectangular box [0,52] x [—40,40] 
with a resolution of 0.8/my, which lead to a grid with 6'666 points. 

The derivatives in the PDEs (jH]) are approximated by finite differences in a straightfor- 
ward way, 

fi^i: ^j) ^ fi,j 5 

{drf){ri, Zj) ^ ifi+ij - fi_ij) , 
{dj){r„ Zj) ^ - fij-i) , 

{dlf){n, Zj) ^ 1 - 2/,, + , 

{drdzj){r i, Zj) ^ (/i+ij+i — — fi-ij+i + fi-ij-i) , 

{dlf){ri, Zj) ^ {fij+i - 2fij + fij^i) , 

where /(rj, Zj) stands for one of the fields Aj., A^p, A^ or (f) at the position (z, j) on the grid. 
Using these replacement rules the PDEs reduce to a system of non-linear equations on the 
internal points of the grid (i.e. i = 2, . . . , — I and j = 2, . . . , — 1). On the axis of 
symmetry r = the finite difference approximation for the PDEs requires special care. The 
regularity conditions (I46cp at r = for the fields Ar and A^ imply 

ArXj = , (Bl) 

= , (B2) 



for j = 2, . . . , Nz — 1. For discretizing Eqs. fl44ap and (144d|) at r = we have to use the 
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replacements 



-dr {t A^p) — IdrA^p 



=0 hi 



1 



=0 hrhz 
2 



{Ar,2,j+1 - ^r,2,i-l) 



=0 hr 



-A 



The system of non-linear equations is completed by the boundary conditions at r = L and 
z = ±Z, which read 



1 and Ai l = Ai^^^ = for z = 1, . . . , A^^ , 



and -^Nrj = for j = 1 



(B3a) 
(B3b) 



Finally, the constraint equation (HSj) is given by 

Nr-l Nz~l 



N, 



cs 



Stt 



EV^ L ^ • T \ A ( ^r,^,j+l - ^r,i,j-l ^z,i+l,j " ^2:,^-l,j ^ 

2_^hr[z-l)A^,,\^ — — J . 



(B4) 

i=2 j=2 

Including the constraint equation ( lB4p we obtain a system of ANrNz + l non-linear equations 
for the unkowns (pij, Aj-^ij, A^^ij, Az^ij and the Lagrangian multiplier fi. This system is 
solved using a standard Newtonian algorithm for non-linear equations. In every step of the 
Newton iteration, one has to solve a system of AN^Nz + 1 linear equations. 
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